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1. Introduction 



The exploration of the phase diagram of QCD and its thermodynamics are challenging prob- 
lems and central goals of lattice simulations at finite temperature and density. See [|l|] for a review. 
In this work we present Monte Carlo results for the thermodynamics of SU (N) gauge theories with 
number of colors, N = 2, 6, in 2 + 1 dimensions. These theories are closely related with those in 
3 + 1 dimensions and are more numerically feasible. We focus on the confined phase T <T C and 
study thermodynamic variables such as the trace of the energy-momentum tensor, pressure, energy 
and entropy density using the integral method. We also investigate scaling properties with N of 
the different observables and compare our results with the predictions obtained assuming that the 
thermodynamics of the system could be described as a gas of free glueballs. We shall show that 
a relevant improvement in the comparison near the critical point is obtained including also higher 
orders in the glueball spectrum and assuming for these terms a bosonic string description. 



2. Thermodynamics on the lattice 

Before discussing the thermodynamics of SU (N) lattice gauge theories in 2 + 1 dimensions, 
we sketch some basic thermodynamics relations in the continuum. From the partition function 
Z(T,V) we get the free energy density as, 

/ = -|io g z(r,y), (2.i) 

where T is the temperature and V is the spatial volume. In the thermodynamic limit the pressure is 
related to the free energy density as, 

p = -\imf. (2.2) 

V-S-oo 

In the following we will assume to have a large, homogeneous system, so that the pressure can be 
identified as minus the free energy. Once the pressure is calculated as a function of the tempera- 
ture p(T), the other thermodynamics variable are derived. For example, the trace of the energy- 
momentum tensor e — 2p is, 

£ ~ 2P -T^-(^]. (2.3) 



r 3 drVr 3 , 

The energy density e = T 2 -^(p/T) is then obtained by adding 2p/T 3 to this result while the 
entropy is given by, 

s= £ -±P = ^. (2.4) 
T dT 

On the lattice the temperature and volume of the thermodynamic system are determined by the 
lattice size N t x N 2 and the lattice spacing a, 

V = (aN s f, T = ±-. (2.5) 

In this work, we perform a non-perturbative study of SU (N) Yang-Mills theories with ,/V = 
2,3,4,5,6 colors regularized on a finite lattice, with lattice spacing a, with N s points along the two 
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space-like directions and N x points along the time-like direction. We use the Wilson action for a 
generic SU (N) gauge group, 

S w (U^x)) = £s(t//>) , S(U P ) =p(l- ^ReTr£/^ , (2.6) 

where P denotes one of the 3N r x N 2 plaquettes on the lattice and Up is the product of the U^- 
matrices (with /I =0,1,2) around each lxl plaquette. On the lattice the partition function is 
given by, 

Z = f n^Mexp(-M^M)) • (2-7) 



In the continuum limit eq. ( J2.6| ) becomes the standard Yang-Mills action provided that, 



t>= 2 4. (2.8, 

In 2 + 1 dimensions g 2 has dimensions of mass and sets the scale. 

Although in principle all thermodynamics variables can be calculated from the free energy 
density, in practice, a direct computation of the partition function on the lattice is not possible. 
Here we use the integral method of Refs. |^, g], as in Ref. [|J]. We first calculate the action, i.e., 
the derivative of the partition function with respect to the bare coupling j8. Up to an integration 
constant, resulting from the lower integration limit j3o, the pressure is then obtained by integrating, 



p(P,N t ,N s ) 



-N 3 T [ P dP'AS(P',N T ,N s ) (2.9) 



T 3 

where in 2+1 dimensions 

AS(P',N T ,N S ) = 3(P ) p - (P s + 2P,) p . (2.10) 

Here P s t denote the expectation values of space-space, space-time plaquettes, respectively and Pq 
is the plaquette value on symmetric lattices A^ 3 . Using eqs. (2.9) and ( 2.10D we can write the trace 
of the energy-momentum tensor (|Q|) as, 



e - 2p - 9 "uW/i( , iw,vc C2.ii) 



In order to obtain eq. ( 2.11 ) as a function of T/T c , where T c is the critical temperature of the 



continuum theory, we need to relate T /T c to j8, for any value of Af, 

p=p(T/T c ). (2.12) 

In practice, T c ^ is determined through a parametric fit (similar to the one performed in [|5|] for 
SU (3)) but for generic A^. A good choice for /3o ^ is j8o = /3 (T /T c = 0.6), after checking from the 
measurements that times the integrand in ( |2. 1 1[ ) is negligible at this temperature. 
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3. SU (N) gauge theories at large N: scaling properties 



Let us investigate the large-V limit of eq. ( |2.12[ ) in SU (N) gauge theories in 2 + 1 dimensions. 
To do so we need to relate some dimensionless ratios that in this limit become constant ^ 0, g|, 
say, 

^=4.108(20) + ^ + ..., ^=0.903(3) + ^+.... (3.1) 



where c is a constant. Here ^/o is the square root of the string tension at zero temperature in the 
continuum theory. 

Considering also that, if we keep g 2 N fixed, j8 scales as N 2 (2JS) and from [§, [7|] 

a 0.12 

0.1975 -—5- + ..., (3.2) 



g 2 N ' N 2 
we get 

r- 0.395V 2 0.24 

y/o = « «-+-■ ( 3 - 3 ) 

ap ap 

Combining these expressions we obtain the dependence of /3 in terms of the temperature T. To get 
j8(r), it is particularly convenient to set the temperature scale using the \fa /T c ratio. To the first 
order in /3 we have 

I = _L^ = r a l (34 ) 

T c T c (0.3957V 2 -0.24) (0.903 + ^) 



and using eq. ( j2.5| ) gives, 



j3 =A^t^ (0.357V 2 + 0.13 -0.2 11 /N 2 ) , (3.5) 



which for N = 3 gives j8 = 0.34 (to be compared with the expression given in Bialas et al. [S], 
J3 = 3.3N r J + 1.5 + 0(1 /V T ), which gives J3 = 0.33). Combining Eq. Q with the data from |g] 
to get the correction to the scaling in the large-TV limit we obtain, 

p =V T ^(0.357V 2 + 0.13 -0.211/V 2 ) + (0.22V 2 -0.5), (3.6) 

which gives the dependence of j8 on the temperature T up to a first order correction to be used in 
eq. (|2~ll|). 



4. Numerical results and discussions 



We are now ready to evaluate the trace energy-momentum tensor in eq. ( [2.1 1| ) and check the 



validity of the scaling dependence in eq. ( p.q ) by plotting the right hand side of eq. ( [2.1 1[ ) vs. 
t = y- This plot is expected not to be dependent on Af (and also on N r ). 

The numerical simulations were performed using the Chroma library JT0| ] plus our own pro- 
grams (for 51/(2) and 51/(4)). We evaluated A5 for V T = 6 (and for SU(2) and SU(3) also for 
V T = 8) and spatial volumes V| such that the aspect ratio was always V S /V T > 8. In agreement 
with Ref. [El], we may safely assume that, in this temperature range, this condition is enough to 
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Trace of the energy-momentum tensor in 2+] dimensions 



Trace of the energy-momentum tensor in 2+1 dimensions 
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Figure 1: Left:The trace of the energy-momentum tensor vs. t = T/T c for SU(N — 2,3,4,5,6). Right: 
Magnified view, of the same, in the low temperature region. 



[race of the energy momentum tensor for SU(3,5.6)vvith SB normalization 



trace of the energy momentum tensor for SU(4,5.6) with SB normalization 
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Figure 2: The trace of the energy-momentum tensor normalized to the lattice Stefan-Boltzmann (SB) limit 
vs. t = T/T c . 



eliminate finite size effects in the spatial directions. A detailed description of our results and of the 



algorithms we used will be reported elsewhere [ |1 3Q 



We report in fig. [I] our estimates for the trace of the energy-momentum tensor. The SU (3) data 
are in perfect agreement with the one in Ref. Below the critical temperature T c there is a good 
scaling with N. Above T c , the different curves split up and they appear to be ordered according to 
the high temperature scaling law for the value of N. See fig. || for the gauge groups SU (3,4,5,6). 

As mentioned in the introduction our main goal was to compare the T <T C data with a glueball 
gas model. We performed this comparison in three steps. First we assumed the gas to be dominated 
by the lowest glueball state, then we included all the glueballs below the two-particle threshold (for 



which very precise numerical estimates exist) and finally, following the suggestion of [11], we try 



to compare our data with the whole glueball spectrum, assuming for the glueballs an ansatz inspired 



by the effective bosonic string model. Details on the calculations can be found in [13]. 

The pressure associated with a single non-interacting, relativistic particle species of mass m 
reads, 

mT 2 1 / . m\ ( T 



2»Ep«n-*fH I+ £;) <41) 



From Eq. (4. 1 ) we can reconstruct all thermodynamics quantities as explained earlier. In particular, 
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the trace of the energy-momentum tensor (2.3) can be written as 



E-2p m 2 y-, 1 / , m\ , A n ^ 

£ -exp (-*- . (4.2) 



r 3 2TiT 2 ^k F V t 

This observable is particularly suited for this comparison since it can be calculated by numer 



ical simulations without being integrated over /3 (see Eq. ( [2.1 1| )). From the above equation we 
obtain: 

m 2 a 2 w-, 1 / m\ 

^-^w^-M-'t)- <43) 

For the first two stages of the comparison we used the numerical values of the glueball masses 
reported in [Q. Given the precision of the data it is mandatory to keep into account scaling correc- 
tions in this comparison. The most effective way to do this is to rewrite m/T as 

- = ^=^- = ^=(aN t V^), (4-4) 

and then use the scaling functions reported in [^]. We write here these corrections explicitly in the 
SU (3) case, for the lowest mass m in the case of a lattice size N T = 8 (the generalization to any 
value of ,/V is straightforward). Using, 

3.367(50) | 4.1(1.7) | 46.5(11.0) ^ (4 f . } 



and -4= = 4.329(41) we obtain, 



m n Z3.367 4.1 46.5 , 

^ = Bx4.329x^ + F + 7r ). (4.6) 



Higher masses can be treated in the same way, using the data for the ratios m,/ yfc reported in [Q]. 
We compare the results of this analysis with the data for the trace of the energy-momentum tensor 
in fig. H for Af = 2 and N >3. The blue and red lines correspond to the inclusion of the lightest 
mass and the first eight masses, respectively. It is easy to see that these fits fail to reproduce the data 
and suggest the necessity of taking into account the full spectrum of glueballs using for instance a 
string inspired ansatz. 

To compare our results with the bosonic string predictions for the Hagedorn spectrum (in the 



same spirit as in Ref . [11] in Ad) we extended to arbitrary dimensions the computation of the 



density of states of the closed bosonic string (following [ 12]). 
We found the following expression, 

7r \^ 1 1 / j ^4-1 ( t h\ m /t h 



~p d -2{M) = [-) -{ d -2)^[^j (4.7) 



Inserting this expression in eq. ( |4.2| ) we found a remarkable agreement with the data, even in the 
region near the critical transition [13]. This comparison is reported in figj3|-left for SU (2) and in 
fig. |-rightforAr>3. 

We think that this type of analysis (which we plan to further improve in the future) will give 
us the opportunity to test the string inspired glueball models (like for instance the Isgur-Paton one 
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Trace of the energy-momentum tensor in 2+1 dimensions 



Trace of the energy-momentum tensor in 2+1 dimensions 





lightest glueball 




first eight 




string prediction 


i 


SLf(2) 





Figure 3: Comparison between contribution of the glueballs spectrum, the string predictions with respect to 
the trace energy-momentum tensor of SU (2) (Left) and SU (N > 3) (Right). 



[]14[]) and also to better understand the many non trivial features of effective string models which 
have been up to now addressed only looking at observables related to the interquark potential or to 
the width of the flux tube. 
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